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Abstract 

We investigate the collective oscillations in the BCS-BEC crossover region of a trapped Fermi 
superfluid at T = 0. We explicitly include the effect of a Feshbach resonance, which leads to 
a tunable pairing interaction adjusted by the threshold energy of the resonance, as well as the 
associated molecules. In a previous paper, we obtained solutions of the Bogoliubov-de Gennes 
coupled equations, describing a trapped Fermi superfluid gas. Using these Bogoliubov quasi- 
particle excitations, we calculate the density correlation function in the Hartree-Fock-Bogoliubov 
generalized random phase approximation (HFB-GRPA). Previous discussions of collective modes 
were based on solving the T = equation of motion for the superfluid in a trap, ignoring the effect 
of single-particle excitations. We determine the frequencies of the quadrupole and the monopole 
modes. In the crossover region, the frequencies of these two modes are shown to be strongly 
depressed by the threshold energy of the two-particle continuum, originating from Andreev single- 
particle bound states (or in-gap states) localized at the edge of the trap. This effect becomes less 
important as one enters the BEC regime, due to the absence of these low energy fermion bound 
states. This suppression of collective mode frequencies by the two-particle continuum may be a 
useful experimental signature of the single-particle energy gap in the crossover region of trapped 
Fermi superfluids. 

PACS numbers: 03.75. Ss, 03.75.Kk, 03.70.+k 
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I. INTRODUCTION 



The recently discovered fermion superfluidity in K[l| and Li 



:Q and 6 Li0,y,y,Q 



has attracted 



much attention. This is the first system where we can experimentally examine the BCS- 



BEC crossover phenomenon 



literature 
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I], [l(3], originally discussed in the superconductivity 
16l.ll7l.ll8||. The novel pairing mechanism of this Fermi superfluid 



in which two Fermi atoms can form a 



is associated with a Feshbach resonance (], [3, _ 
quasi-molecular (boson) state. The resulting attractive interaction can be tuned by varying 
the threshold energy (denoted in this paper by 2u) of the resonance, using an external 
magnetic field (21L l^. 

In Fermi superfluids, there are two distinct kinds of excitations, i.e., collective modes 
and single-particle excitations. Two single-particle excitations are created with the breakup 
of a Cooper-pair. Even in the BEC regime, where the character of superfluidity can be 
best described as the Bose condensate of molecules, the single-particle Fermi excitations 
still exist. This situation quite contrasts to an ordinary BEC of Bose atoms, where low- 
energy excitations are all bosonic in nature 2^ . Very recently, both the single-particle 



excitations 



2J] and collective modesjj, |25( have been studied experimentally in the crossover 



region of superfluid 6 Li. 

The frequency of a collective mode is strongly affected when it is near the threshold energy 
of the two-particle continuum at 2E g , where E g is the single-particle excitation gap. In the 
BCS region of a uniform, neutral Fermi superfluid, this can be seen by the suppression 
of the Goldstone mode 1^ = ^9 at low q), when its energy approaches the two-particle 
continuum value of 2|A| 2f|. Here, Dp is the Fermi velocity and A is the BCS Cooper-pair 
order parameter. A similar suppression of a collective mode frequency by the single-particle 
excitations is expected in a trapped Fermi superfluid. In this case, we have shown recently 2]| 
that the single-particle energy gap E g is not simply related to the position-dependent order 
parameter A(r), but is dominated by the lowest energy Andreev bound state appearing at 
ihe eogeot the I trapped gaS Q. Even in the crossover region, where the magnitude of 
the superfluid order parameter in the center of the trap is comparable to the Fermi energy 
e-p, Eg is still the order of the trap frequency ujq (which is usually much smaller than £p) . 
Since the frequencies of the quadrupole and the monopole modes are the order of cl>o[30], 
they are expected to be strongly affected by the two-particle continuum at 2E g . 
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In this paper, we investigate the collective modes at T = in the BCS-BEC crossover 
region of a Fermi superfluid in a trap. The theory is precisely the same as in our recent 
discussion of collective modes in a uniform Fermi superfluid^], except that we limit ourselves 
to T = 0. The present paper makes crucial use of our results [27] for the single-particle 
excitations in a trapped Fermi gas superfluid. To properly include the effect of the energy 
gap Eg associated with the Andreev bound states, the self-consistent Bogoliubov de Gennes 
(BdG) coupled equations 3l| were solved numerically. The strong-coupling effect is taken 
into account by considering the deviation of the chemical potential \x from the Fermi energy 
£fj which is associated with the appearance of the bosonic pair states. As a result, the BdG 
equations have to be solved together with the equation of state involving the number of 
atoms [27] to take into account the molecule condensation. 

In this paper, the single-particle solutions of the BdG equations are used in a detailed 
calculation of the density correlation function, using the time-dependent Hartree-Fock- 
Bogoliubov (random phase) approximation (HFB-GRPA). This includes the coupling be- 
tween ordinary density fluctuations and fluctuations in the Cooper pair-channel. As proven 



theorem 



recently 32 , this linear response formalism guarantees that Goldstone's theorem and Kohn's 



33] are satisfied in the entire BCS-BEC crossover regime, required for any consistent 



approximation for collective modes in a Fermi superfluid. 



3J, |35| to the 



Our approach is an extension of the weak-coupling BCS theory by Bruun 
strong-coupling regime including a Feshbach resonance and the associated molecules. In 
the crossover region, we show that the frequencies of the low-energy collective excitations, 
such as the quadrupole and monopole modes, are strongly suppressed so that they stay 
below the threshold energy of the two-particle continuum at 2E g . We note that the T = 



hydrodynamic superfluid theory 



3(1 36, 



37 



381 ] describing pure superfluid dynamics cannot, 



in principle, describe this suppression phenomenon. This is because such a theory ignores 
any explicit consideration of single-particle excitations, and in particular the low energy 
Andreev excitations in a trap. As expected, this suppression effect disappears as one enters 
the BEC region, where the single-particle excitation energy gap becomes very large [271] . 

It is useful to be clear about the relation of our kind of microscopic theory of collective 
modes and the predictions based on solving the equation of motion of a pure superfluid 3(3, 



36 



38j . The latter so-called "hydrodynamic" theory is limited to T = 0, since it assumes that 



all the atoms are described by the superfluid component (i.e., there is no normal fluid). In 
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this theory, the only unknown microscopic input is the equilibrium equation of state, giving 
the chemical potential as a function of the density. There has been considerable recent 
progress in determining this equation of state in the BCS-BEC crossover regime, including 



39, 



4Cj . The advantage of starting from 



the strong interaction unitarity limit ((/cpap) - = 0) 
the equation of motion for a pure superffuid is that it is expected to be valid even in the 
strong interaction limit, and makes no assumption about the amount of depletion of the 
condensate (which could be very large in the unitarity limit). The weakness of the T = 
superfluid hydrodynamic theory is that, by definition, it cannot deal with any questions 
related to the single-particle excitations, such as the effect of the two-particle continuum 
or the effect of a normal fluid component (which would arise at finite temperatures). In 
contrast, our present microscopic calculations can deal with such effects. However, they do 
not include any large depletion of the condensate due to interaction, although the effect of 
thermal excitations can be treated. 

One deficiency of our present calculations might be mentioned at the outset, namely 
that in the BEC limit, our collective modes are predicted to be those of a non- interacting 
molecular Bose condensate. The explanation of this is simply that in mean-field BCS the- 
ory which lies at the core of our present discussion, we do not include the higher order 
fluctuations which ultimately give rise to the correct s-wave interaction between molecules 



421 ] . we give the 



cim = 0.6ap|41(, in the BEC limit of the crossover. In a separate paper 
appropriate extension which does reproduce the correct value of ayi, which in turn would 
lead to correct collective mode frequencies in the BEC limit j3,Q|- However, this deficiency 
in reproducing the BEC collective modes correctly is of minor concern in connection with 
the major point we emphasize in this paper, namely the suppression of the collective mode 
frequencies in the BCS-BEC crossover region due to coupling to the two-particle continuum. 

It should also be emphasized that our microscopic HFB-GRPA theory for the collective 
modes could be used as the starting point to derive the equation of motion for the the 
dynamics of the superfluid at T = 0. In the BCS limit of a trapped Fermi superfluid, for 
example, such a derivation has been carried out in detail by Baranov and Petrovj^j]. The 
analytic solutions of this T = equation of motion describ ing the collective modes of a BCS 
Fermi superfluid have been discussed by Bruun and Clark[46|. 

This paper is organized as follows: In Sec. HH we briefly review^, 27] our coupled 



fermion-boson (CFB) model in the presence of an isotropic harmonic trap. We calculate 
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the generalized density correlation functions based on the HFB-GRPA formalism in Sec. 
TTT1 In Sec. IIV1 we summarize our previous results[27] for the single-particle excitation 
spectrum in a trapped superfluid Fermi gas, and then show in Sec. Mhow the single-particle 



energy gap depresses t 



mentioned in Refs. 



27 



re frequency of collective modes. Some of our results were briefly 



471 ] . Appendix A gives a brief review of the results for the single- 



particle BCS excitations obtained in Ref. |27j. These are the input for the calculation of 
response functions in this paper. Appendix B proves that the renormalized Bose propagator 
for bound states has the same poles as the fermion density response function, as expected. 



II. COUPLED-FERMION BOSON MODEL 

We consider a two-component superfluid Fermi gas with a Feshbach resonance trapped in 
a harmonicpotential. To describe this, we consider a coupled fermion-boson (CFB) model 



given by 



mc potential. 

flflflflfl 



27 



V 2 



H = Y,j »t(r)[-— + ^ ap (r)-//]* CT (r)-C/y dr*f(r)*I(r)* ; (r)* T (r) 



+ / dr$ f (r 



2M 



+ 2v + CW - 2 ^l $ ( r ) + 9r / dr[$t( r )* ; (r)* T (r) + h.c. 



(2.1) 



This model involves two field operators for Fermi atoms ^/ (T (r) with pseudo-spin a =1,1 
and a Bose field operator $(r) associated with the Feshbach resonance. The Feshbach 
resonance with strength g r describes the conversion of two Fermi atoms into a molecular 
boson and the dissociation of a molecule into two Fermi atoms. This resonance is controlled 
by the threshold energy 2v of molecular excitations, which is assumed to be adjustable. U 
represents a non-resonant interaction, which we take to be attractive (-U < 0). Since a 
molecule consists of two Fermi atoms, we take M = 2m and impose the conservation of the 
total number of atoms, 

N = N F + 2N M = ]T / dr<tft (r)* ff (r)> + 2 f dr($t(r)$(r)>. (2.2) 

This constraint has already been taken into account in Eq. (J2.1|) by assuming that Bose 
chemical potential \im = 

where \i is the Fermi chemical potential. In this paper, 
we assume that the harmonic trap potentials for Fermi atoms V£ (r) and Bose molecules 



^tra P (r) are both isotropic, as 



trap 



Cap = \Mul M r\ 



(2.3) 



2 u > trap 2 

We also assume that atoms and molecules feel the same trap frequency (uom — ^>o), which 
describes current experiments using Feshbach resonances. 

Since the CFB model in Eq. (|2.1|) involves Fermi atoms and Bose molecules, superfluidity 
is characterized by the BCS Cooper-pair potential A(r) = (^j(r)\I/^(r)) and the molecular 
BEC condensate 4>M (r) = ($(r)). In equilibrium state, these order parameters are related 
to each other by|27j] 

V 2 



9r_ 

U 



A(r) + 



2M 



+ Or)+2i/ -2/1 



0. 



(2.4) 



This identity ensures that in the superfluid phase, A(r) and <pM{r) are both finite in the 
same spatial region and vanish together. 

In discussing collective modes, it is convenient to divide Eq. (j2.1|) into the Hartree-Fock- 
Bogoliubov (HFB) mean-field Hamiltonian (= -£^hfb) and a fluctuation part (= Hp^)[26\. 



We introduce the Nambu representation for fermion field 

\ 



49] 



tt T (r) 
*t( r ) } 



(2.5) 



and its conjugate field ^(r) = (\l/j(r),\l/ 
mean-field Hamiltonian -£/hfb reduces to[27| 



r)). In this Nambu representation, the Fermi 



HFB 



U , s 
— - -n F (r) + V t 



trap V 



r 3 - A(r)n 



*(r). (2.6) 



2m 2 

Here Tj (i = 1,2,3) are the Pauli matrices acting in particle-hole space. np(r) 



X] (T (^ / j T ( r ) l ^o-( r )) i s ^e number density of Fermi atoms. We note that the fermion part 
i^HFB nas t ne same form as the standard BCS Hamiltonian except that the Cooper-pair 
order parameter A(r) is replaced with the composite order parameter, given by 



A(r) = A(r) - g T (t> M {r) 



(2.7) 



In Eq. (|2.6|) . the order parameter A(r) is taken to be real and proportional to the T\- 
component, without loss of generality. Because of the assumed spherical symmetry of the 
trap, A(r), A(r), 4>m ( r ), and n-p{r) are all isotropic depending only on r = |r|. 



The molecular Bose Hamiltonian H M is given by the third term in Eq. (|2.1jl . The part 

5H M = J dv8^(v)[-^ + 2u + V^ p (r)-2p]S^(v) (2.8) 

describes fluctuations from the Bose condensate, where the molecular quantum field operator 
is defined as <5$(r) = $(r) — (f) M {r). 

Fluctuation effects around the mean-field theory described by -£^hfb can be expressed as 
interactions between fluctuations in the density-channel (particle-hole channel) and in the 



Cooper-channel (particle-particle channel) [26] . It is convenient to introduce the generalized 
density operators, defined by 5(| 

Pj (r) = &{t)tM*) U = 1.2,3). (2.9) 

Since we take A(r) proportional to the ri-component, pi(r) and p2(r) represent amplitude 
and phase fluctuations in the Cooper-channel, respectively. p 3 (r) is the ordinary density 
operator, expressed in the Nambu representation. Using Pj(r), we can write the non-resonant 



interaction U as 



5l| 



-- J dvpj^pjiv) (7 = 1,2,3). (2.10) 

Here Hp^ and HpQ involve the the Cooper-channel, while H F 3 ^ involves the density- channel. 
The Feshbach resonance or the last term in Eq. ()2.1|) can be expressed as an interaction 
between fluctuations in the molecular field <5<3>(r) and the Cooper-pair fluctuations pi(r) and 
p 2 (r), namely 

= 9 - fdr 



J dr[5$ f (r)p_(r) + <5$(r)p+(r)] , (2.11) 



2 

where we define p±(r) = pi(r) ± zp2(r). Eq. ()2.11|) shows that the molecular fluctuations 
5$(r) couple with the fluctuations pi(r) and P2(r) in the Cooper-channel, but do not couple 
with ordinary density fluctuations described by ps(r). 

The fermion mean-field Hamiltonian H^ FB in Eq. (|2.fijl can be diagonalized, redu cing 
the problem to solving the Bogoliubov-de Gennes (BdG) coupled equations. In Ref. |27| . 
we have self-consistently solved these BdG equations and used the results to calculate the 
static properties as well as single-particle excitation spectrum in the BCS-BEC crossover, 
at T = 0. Since the single-particle excitations are the ingredients in our calculations of the 
density correlation function, we summarize some of the discussion about BdG equations in 
Appendix A. For more details, see also the discussion in Ref. j^J. 
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III. GENERALIZED DENSITY CORRELATION FUNCTIONS 



Collective modes corresponds to the poles in the density correlation function. In 
the superfluid phase, density fluctuations couple with particle-particle or Cooper-pair 
fluctuations 26], so that one needs to take their coupling into account. For this purpose, it 
is convenient to introduce the generalized density correlation functions, defined by 

U ij {ry,u) = -i / dte^Hpiir^Jjir'^)}). (3.1) 

J 

n 33 is the usual density correlation function. Iln and J\ 22 are the amplitude and phase 
correlation functions, respectively. The off-diagonal components 11^- describe couplings 



between different kinds of fluctuations. For examp 



to the Josephson effect, as discussed in Refs. |8 



e, n 2 3 is the phase-density coupling due 



26] 



We calculate the generalized density correlation functions Ily in a trap, based on the 
HFB-GRPA formalism in terms of the interactions Hp^ in Eq. ()2.10|) and in Eq. 

()2.11|) . This is simply an extension of our previous work for a uniform superfluid Fermi 
gasjs| to include an isotropic harmonic trap. In contrast to Ref. jsj], this paper considers the 
case T = o[^. We refer to this earlier paper for a more detailed discussion of the physics. 
The 3x3 matrix correlation function fl = {Ily} (i,j = 1, 2, 3) defined in Eq. (|3.1|) satisfies 
the GRPA equation in real space, namely 

fl(r,r» = fl°(r,r» - ^ / dr"fl(r, r", u)fl°(r", r, u) 

+ 9 1 J dr" J dr'"n(r, r", u)B(r", v'" , u)Ti\r", r', uo). (3.2) 

In a matrix notation for the dependence on position, the formal solution of Eq. ()3.2|) is 
given by 

'"tf-^sMifivyr 1 . (3.3) 



1+ 2 



In Eq. (13.2)1 . the effects of a Feshbach resonance appear as the non-local and energy- 
dependent effective interaction B(r", r'", uj). 

B has only the matrix elements B u , B 12 , B 2 i, and B 22 , indicating that it does not 
involve the ordinary density P3(r). This originates from the fact that the Feshbach resonance 
is described by an interaction between fluctuations of the molecular Bose field <5$(r) and 
Cooper-pair order parameter fluctuations [see Eq. (|2.11j) ]. These non-zero matrix elements 
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can be expressed as 



B n (r,r',u) B 12 (r,r',u) 
B 2 i(r,r',uj) B 2 2(r,r',u) 



PD°(y,y\uj)P\ 



(3.4) 



where the 2x2 matrix P is defined as 



V2 



1 1 

i —i 



(3.5) 



and £>°(r,r',w) is the 2 x 2-matrix molecular Bose Green's function given by 

fM l \ fM ( /\ 
J nlm\ J J nlm\ ) 

u + t 3 - (Eg + 2u- 2/i) 



nlm 



l in 



where the l-th component has the explicit form, 



(3.6) 



(3.7) 



u + T 3 -(E% + 2v-2ti)' 

In Eq. (|3.2|) . II (r,r / ,u;) = {n? } is the zeroth order 3x3 matrix density correlation 
function, obtained from the analytic continuation of the corresponding two-particle thermal 
Green's function, 



1 



n°-(r, r', iv m ) = - J2 Tr [TiG(T, r', iv m + iu; m )r i G(r / , r, iu m ) 

P Urn 



(3.8) 



Here, iu m and iv m are the fermion and the boson Matsubara frequencies, respectively. The 
2x2 matrix single-particle Green's functions have the form, 

G(r,r',^ m ) = - jfVe^ T (r T {tf(r,r)*V>0)}). (3.9) 

Using the eigenfunctions fi mr Sj) of the harmonic potential, we can write Eq. (J3.9J) 

G(r,r',^ m ) = J^Yirnim^r^i^Y^e'), (3.10) 

lm 

in terms of the fermion Green's function g l (r, r', iuJ m ) for angular momentum I, 



g l (r, r',iu m ) 



E 



3 L 



Mr)A^(r') A 3 - t (r)A;,(r 



t / 



+ 



Here we have introduced the wavefunctions related to BdG solutions 

W l 

W N t +n,j+l 



(3.11) 



(3.12) 



Substituting Eq. (|3~TUJ) into Eq. (J31|, one obtains Q 

n°(r,r> m ) = £Y iM (0>°/(r,r> m )Y L * 



(3.13) 



LM 



Here, n° L (r, r', it/ m ) describes the (generalized) density fluctuations with angular momentum 



n°f(r,r'> m ) 



(3.14) 



„, 47r(2L + l) 

where (101'0\L0) are the Clebsch-Gordan coefficients [53]. This leads to a selection rule in 
the summations with respect to I and I'. The r-dependence of the correlation function is 
determined by rfJr,r' ,iu m ), which has the form 



P Urn 



Ti§ l (r, r', %v m + iu} m )Tjg l (r', r, iu r< 



= ^Tr[r,A ni (r)AL(rOr,A n ,(r')Ai,(r) 

nn' 



1 



1 



iv m + (E% + E&) ' 

The correlation function r"' (i = 1, 2, 3), and T 1 ^ can all be written in the form 



nn' 
M 1 



rg (r, r', w ro ) 



+ (E nl + E nn 



F \2' 



where 



Q""'(r,r') = Trkr,A ni (r)AL(r')r J r 2 A n Kr')AL(r; 



The remaining components are given by 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



After carrying out the usual analytic continuation, iz/ m — > u; + i8, in Eq. (|3.13p . we 
substitute Eqs. (|3.6|) and ()3.13|) into Eq. (|3.2)l . The result for Il(r, r', uj) can be decomposed 

as 

(3.19) 



J2Y LM (6)U L (r,r',uj + )Yl M (e'), 

LM 



where the L-th component Ii L satisfies the GRPA equation 

U L (r, r' uj) = U 0L (r, r' uj) - - r" 2 dr"U L (r, r", uj)U 0L (r'\ r' uj) 



rp' POO pOG . 

+ — / r" 2 dr" r"' 2 dr'"n L (r,r",uj)B L (r",r , ",uj)n 0L (r",r , ,uj). 
2 Jo Jo 

(3.20) 
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B L (r, r', u) is given by Eq. (J3.4j) but now with the Bose propagator D°(r, r', to) replaced with 
its L-th component D 0L (r,r' ,u) defined in Eq. ()3.7|) . In matrix notation for the position 
dependence, the formal solution of Eq. ()3.20|) is given by 

1 r 



n 



OL, 



OJ 



1 + 



U-g 2 B L (uj 



it 



OL, 



OJ] 



(3.21) 



where U L (r, r',u) = rU L (r,r' ,u)r' , U 0L (r,r' ,u) = rtl 0L (r, r', oj)r', and B L (r,r' ,u) = 
rB(r, r', oj)r' . 

The frequency of the collective mode in the L-th channel can be obtained from the pole 
of the response function U L (r, r', u). The equation for the energy of the collective mode is 
given by the zero of the determinant 

1 



det 



1 + 



U-g£B L (u 



IT 



OL , 



0. 



(3.22) 



The dipole, quadrupole, and monopole modes are described by the L — 1, L — 2, and L = 
channels, respectively. All the collective modes accompanied by density fluctuations and/or 
fluctuations in the Cooper or particle-particle channel are given by the solutions of Eq. 
(I3.22J) . For example, the gapless Goldstone mode at uo = appears in the L = channel. 
We note that molecular boson Green's function has the same poles as those obtained from 
the solutions of Eq. ()3.22|) . as proven in Appendix B. The physics behind this equivalence is 
that in the coupled fermion-boson model, fluctuations in the molecular field strongly couple 
into the fluctuations in the Cooper-channel through the Feshbach resonance. Thus, we need 
only calculate the density correlation function for fermions to find the collective modes. 

The frequencies of the lowest energy collective modes are determined numerically by 
solving Eq. (j3.22j) . The results will be discussed in Sec. El A weighted strength function 
for various density fluctuations is usefully defined as 



S 33 (u,L) = / F L {r)r 2 dr / F L {r')r' 2 dr'Im U^ 3 (r,r',u) 

n Jo Jo L 



(3.23) 



where Fi(r) = r for the dipole mode, ^(r) = r 2 for the quadrupole mode, and Fo(r) = r 2 
for the monopole mode[54]. These strength functions are evaluated in Sec. El 
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IV. SINGLE PARTICLE EXCITATIONS IN A TRAPPED SUPERFLUID FERMI 
GAS 



In a recent paper [27j, we have numerically solved the BdG coupled equations together 
with Eq. ([A 13)1 . which determines the Fermi atom chemical potential /i. We self-consistently 
determined the Cooper-pair order parameter A(r), the molecular condensate <$m {f), the 
atom density n F (r), and the chemical potential fi. In the crossover region, we showed 
that the single-particle excitation gap E g ~ uq is much smaller than the magnitude of the 
composite order parameter A(r = 0) ~ ep ^> ojq at the center of the trap. Since this is the 
key point in understanding the behavior of the collective mode frequencies in the BCS-BEC 
crossover, we briefly summarize the single-particle prop erties in a trapped Fermi superfluid 
in this section. The details are worked out in Ref. [27J and reviewed in Appendix A. 

In our numerical calculations, we take N = 10, 912, which is equivalent to filling up states 
up to the energy E = 31.5fru;o (= £f) for a non-interacting Fermi gas in a trap. The Feshbach 



coupling is taken as g r /yRp = 0.06cj (#r = 9ryN/Rp = 0.2e F , where Rp = \j2E-p/mujQ 
is the Thomas- Fermi radius for a free Fermi gas in a trap). In order to see the effect of 
the single-particle excitation gap E g on collective modes, we consider the two cases for the 
nonresonant interaction U: U/Rp = 0.001a;o (this corresponds to U = UN/R F = 0.35£f) 
and 0.0015a;o (U = 0.52ep)- In the BCS region, the choice U = 0.35ef gives E g ~ ojq and the 
case U = 0.52£ F gives E g ^> tj n , where u is the trap frequency. For the high energy cutoff, 



we take tu c = 161.5^0 (3> £f)|55J- in order for the theory to be cutoff-free, we introduce 
the scattering length a s defined by[?], 27] = U c s/[1 — C4ffEo,o; c 2i~]> where e p = p 2 /2m 
and U e s = U + g 2 /(2u + (3/2) ^a;o — 2/x) is the effective pairing interaction. For further 
discussion, see Ref. |27j ]. 

The single-particle density of states N(u) is defined by 



N(u = / drlm G u (r, r, iu m — > 10 + iS) , (4-1) 

71 J L J 



where Gu is the (l,l)-component of the (analytic-continued) matrix single-particle Green's 
function defined in Eq. (|3.9J) . From the spectrum shown in Fig. 1, one sees that a finite 
excitation energy gap occurs at E g ~ 1.2u;o. The threshold energy to break a Cooper-pair 
is given by 2E„. In a uniform BCS superfluid, the Bogoliubov quasi-particle excitation 
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spectrum is given by the well-known expression 



S p = V( £ p-^) 2 + I A I 2 - (4-2) 

In the weak-coupling BCS regime, the chemical potential is close to the Fermi energy of 
a non-interacting Fermi gas, \i ~ e-p > 0. In this regime, the energy gap E g is equal 
to the order para meter, E g = |A|. This relation is no longer valid in the strong-coupling 



ord er p ara 



regime 112. [I3j,|l5|. When the chemical potential becomes negative, we find E g = \J jj? + |A| 2 . 
Moreover, the relation E g = |A| is not valid even in the BCS regime when we deal with a 
trapped gas. Indeed, as shown in Fig. |21 the position-dependent composite order parameter 
A(r) = A(r) — g r <t>M{r) is nearly equal to lOt^o at the center of the trap, which is much 
larger than the value E g = 1.2c<Jo shown in Fig. ^ Even at r = 0.5-Rf, the order parameter 
is still much larger than E g , with A(r = 0.5i?p) ~ 5u 3> E g . 

The difference between E g and A(r = 0) is most remarkable in the crossover region. As 
shown in Fig. E^a), the value of A(r = 0) simply increases as one increases the strength 
of the (Feshbach-induced) pairing interaction. However, the energy gap E g is insensitive 
to this increase of the pairing interaction, always being of the order of the trap frequency 
ujq. When the chemical potential becomes negative in the BEC regime [see panel (b)], E g 
increases, approaching E g ~ |^||27|. In this BEC regime, 2E g physically gives the threshold 
dissociation energy of a tightly bound molecule, i.e., 2E g = \(j,m\ = 2|/i|. 

The reason why the values of E g and A(r = 0) are very different is because of the 



Andreev bound states appearing at the edge of the trap 



In a trapped superfluid, atoms 



feel the "off-diagonal" pair potential A(r) (solid line in Fig. EJ) as well as the "diagonal" 



trap potential measured from the chemical potential 56] 



VS*(r) = CV£p(r) - A*)e(Cp(r) - //), (4.3) 

(dashed line in Fig. |2J), where Q(x) is the step function. Although A(r) and V^g ap (r) appear 
in the diagonal and off-diagonal components, respectively, in the BdG coupled equations, 
they can be thought of roughly as giving rise to a combined potential well 

U^A(r)+CW. (4.4) 

As a result, the Andreev bound states (or in-gap states) appear well below A(r = 0), being 
localized around the bottom of this combined potential well. The lowest bound state has 
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energy E g . This phenomenon is analogous to the appearance of surface bound states in 
superconductivity, originally discussed by de Gennes and Saint James 2t| • There is a large 



literature on t 
As in Ref. 



rese states. 



2jj, it is useful to also calculate the local density of states N(u, r), defined by 

N(u, r) = Im Gn(r, r, iu m — > uj + iS) . (4.5) 

Figure 4(a) exhibits the discrete in-gap (or Andreev) states at frequencies well below the 
order parameter at the center, uj < A(r = 0) ~ lOcuo- These low energy states are localized 
at the edge of the trap r ~ Rp. Such localized surface states are not given by a local 
density approximation (LDA) [see Fig. EJb)]. In the LDA, the local density of states 
N(u, r) is simply given by the uniform BCS density of states, with the order parameter 
and the chemical potential being replaced by the position-dependent ones, A(r) and fi(r) 
// — V£ (r), respectively. The ridge peak in the LDA density of states in Fig. 4(b) simply 
corresponds to A(r). 

We note that the position of an in-gap state depends on the magnitude of the off-diagonal 
pair potential A(r) at the center of the trap. We have seen that the bound state wavefunction 
is well localized at the edge of the trap when A(r = 0) ^> ujq. On the other hand, the surface 
bound states penetrate deep inside the trap when A(r = 0) is small. In Fig. E[a) [where 
A(r = 0) ~ 1.86u;o], the local density of states at the center of the trap (r = 0) has a large 
intensity below uj = A(r = 0) ~ 1.86u;o, originating from the Andreev bound state. 

In contrast, in the BEC regime [(/cpa s ) _1 < 0], where the chemical potential \x is negative, 
the low-lying bound states at the edge of the trap become less important for the single- 
particle excitation spectrum. As expected from Eq. ()4.2j) for a uniform Fermi gas, the 
threshold of the excitation spectrum is then dominated by the chemical potential \\x\. There 
are no Andreev bound states in the local density of states shown in Fig. E(b). Instead, a 
very large excitation gap appears, due to the large dissociation energy (~ 2|/x|) of a molecule. 
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V. COLLECTIVE MODE FREQUENCIES IN A TRAPPED FERMI SUPER- 
FLUID 



A. Quadrupole mode 



In the T = hydrodynamic theory appropriate to a pure superfluid, a collective mode 
of a superfluid Fermi eras in the BCS weak-coupling limit in an isotropic harmonic trap has 



the frequency 



uu = u Jl+ -n(2 + L + n), (5.1) 



where L and n are the angular momentum and radial quantum numbers of the collective 
mode, respectively. This expression predicts that the frequency of the quadrupole mode 
(L = 2, n = 0) is lo = \/2ujq. Figure EH shows the frequency of the quadrupole mode in the 
BCS-BEC crossover, as given by our calculations summarized in Sec. IIHI The quadrupole 
mode frequency is suppressed so as to stay below the two-particle continuum threshold 
given by 2E g , where E g is the single-particle energy gap calculated in Sec. IIHI In the BCS 
regime, our numerical results in Fig. 6(a) for the quadrupole mode agree with the prediction 
to = \/2ujq. In this case (U = 0.52£f), the threshold energy of two-particle continuum 2E g 
is greater than the collective mode frequency y/2io . However, when 2E g < V2coo, Fig. Efb) 
shows a strong deviation from the predicted hydrodynamic frequency to = \[2loq even in the 
BCS regime. We note that quadrupole mode in our calculation has no damping (due to 
decay into two quasiparticles) in the entire BCS-BEC crossover region, since the frequency 
always lies below the computed values of 2E g . 

The suppression of the quadrupole mode frequency due to the two-particle continuum 
shown in Fig. HJb) continues until we enter the BEC regime, where the chemical potential 
becomes large and negative and hence 2E g also becomes large. In Ref. j^f, we showed that 
in the BCS-BEC crossover region, the combined potential well in Eq. ()4.4|) becomes wider 
in the crossover region. This leads to a decrease in the energy of the lowest Andreev level, 
which in turn determines the single-particle energy gap E g . Reflecting this phenomenon, 
in Fig. E[b), we note that the suppression of the quadrupole mode frequency is strongest 
around (kpa s )~ l ~ 0.5, where the value of E g is a minimum. 

Since the single-particle excitation gap E g becomes very large (^> co ) when /z < (see 
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Fig. |3J), the effect of the two-particle continuum on the quadrupole mode is weak in the 
BEC regime. As a result, in Fig. E^b), the frequency of the quadrupole mode shows a sharp 
increase as we enter this region. In the extreme BEC limit [(/ci?a s ) _1 — > oo], our quadrupole 
mode frequency approaches the energy level of a non-interacting molecular Bose gas with 
L = 2 and n = 0, which is given by uj = (2n + L)uj = 2uj . Thus our present calculations 
do not give the correct quadrupole mode frequency in the BEC limit for an interacting 
Bose-condensed gas of molecules j^j]. We briefly explain the reason for this. 

In our present approximation, the effective interaction between molecules (= Vjg ) in the 
BEC regime is mediated by free Fermi atoms, as shown in Fig. In a uniform superfluid 
Fermi gas, taking incident and scattered molecules to be at zero momentum, this process 
gives the following effective interaction 

V# = f E Gl x {v,iu m )Gl x {-v,-iu m ) = ^!|t^W (5-2) 

It is interesting to note that this expression for the repulsive molecular interaction is precisely 
in agreement with that computed from the velocity of the Goldstone phonon in the BEC 
region of a uniform Fermi superfluid (see Eq. (4.20) of Ref. Q). As we have discussed 
earlier, in a trapped Fermi superfluid, our present calculations are limited to a weak Feshbach 
resonance, where the coupling parameter g x = 0.2sf < £p, due to computational problems. 
In this case, by analogy to the uniform gas result in Eq. (|5.2jl . we expect Vjg (oc gf) to 
be very small. This explains why in our results in Fig. |HJ the quadrupole mode frequency 
approaches the non-interacting Bose gas result uj = 2ujq in the BEC region dominated by 
molecules (where 2E g is large and \i < 0). 

For a broad Feshbach resonance (g T ^> ep), in contrast, one can expect a stronger effective 
interaction Vj^ between molecules. In this case, in the BEC regime, we would exp ect that 



the collective modes would be those of trapped Bose superfluid at T = 



30 



|43j, namely 



uo = luoVL + 3n + 2nL + 2n 2 . (5.3) 

This T = hydrodynamic superfluid theory gives uj = \/2ujq for the quadrupole mode in 
the molecular BEC region, the same as in the BCS limit (in contrast to a uniform gas). 
Since the threshold energy 2E a becomes small near the superfluid transition temperature 

n d 

T c [6, 7jJ, the suppression of the quadrupole mode frequency by the two-particle continuum 
will become more important at finite temperatures. Even in the case shown in Fig. Efa), 
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the quadrupole mode frequency will deviate from the simple prediction to = \/2uj in the 
BCS and the crossover region, as soon as the threshold energy of the two-particle continuum 
2E g (T) is greater than V2u . Thus, observing the suppression effect on the quadrupole 
mode frequency by the two-particle continuum might give useful information about the 
temperature dependent single-particle energy gap E g (T). 

Figures |H1 is a plot of the strength function £33(0;, L = 2) defined in Eq. (|3.23J) for the 
parameters used in Fig. 6(b). In each panel, the lowest energy peak is the quadrupole mode. 
In panels (a)-(c), since the threshold energy of the two-particle continuum (shown by the 
arrows) is close to the frequency of the quadrupole mode, we can also see fine structure 
at frequencies just above the quadrupole mode, associated with the two-particle excitation 
spectrum. On the other hand, in the case U = 0.52ep, the threshold energy 2E g of the two- 
particle excitation spectrum is much higher than the frequency of the quadrupole mode, as 
shown in Fig. 6(a). As a result, only a peak associated with the quadrupole mode appears 
in the low-energy spectrum of the strength function, as shown in Fig. 9. There is no fine 
structure associated with the two-particle continuum, of the kind shown in Fig. |H1 



B. Monopole mode 



Figure 10 shows the frequency of the monopole mode in the BCS-BEC crossover region. 
This collective mode is also suppressed below the two-particle continuum at 2E g . (Note 
that the BCS superfluid hydrodynamic prediction in Eq. (|5.1j) gives uj = 2uj for L = and 
n = 1). Thus, as in case of the quadrupole mode, the monopole mode is not damped at 
T = 0. In the BEC regime (fj, < 0), the effect of two-particle continuum on the monopole 
mode frequency is absent because 2E g ^> u . In the BEC limit, our calculated monopole 
mode frequency approaches to = 2u)o, the excitation frequency for a non- interacting Bose 
gas with L = and n — 1. The quadrupole mode has the identical feature, for reasons 
discussed above. When the effective interaction between molecules is significant, as expected 
in the case of a broad Feshbach resonance, the monopole mode frequency would approach 



u = y/Eooo — 2.2Aujq in an improved theory |4S 

The monopole mode appears as the lowest energy peak in the spectrum of the strength 
function Szz{u,L = 0) for the L = channel, as shown in Figs. ITT1 and fT2l As expected, 
we can see the spectral weight coming from the two-particle excitations near the monopole 
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peak shown in panels (a)-(c) of Fig. ^2 

The situation is quite different in the dipole (Kohn) mode|32|,|33(. This collective mode is 
the rigid oscillation of the center of motion, so that the interaction between atoms, dependent 
only on the relative coordinate of two atoms, has no effect on the frequency of this mode. 
As a result, the dipole mode always appears at the trap frequency u , irrespective of the 
value of the two-particle continuum 2E g . 

C. Comparison to experiment 

In concluding this section, we would like to briefly discuss the possible relevance of our 
prediction to current experimental results on collective modes in trapped Fermi superflu- 
ids. The radial breat hing modes have been studied experimentally in a 6 Li superfluid in 
a cigar shaped trapjj, 13] . The observed modes have small damping (r) in the BEC and 
crossover reg ion, but T suddenly increased precisely where frequency of the mode suddenly 
increased j2a|. As the unitarity limit [(pp^s) -1 = 0] is approached from the BEC side, 
the radial breathing mode frequency always decreases, in contrast to the prediction by the 
T = superfluid hydrodynamic theory ffifl. At the unitarity limit, the frequency of the 
radial breathing mode in Ref. [l| agrees with the "universal" value unitarity) predicted 
by superfluid hydrodynamic theory. However, in Ref. j^j], the observed radial breathing 
mode frequency was found to be smaller than ojunitarity Although our calculations are for 
an isotropic trap, the suppression effect we obtain for the collective mode frequencies by the 
two-particle excitation threshold should also arise in anisotropic traps. This is ignored in 
the T = equation of motion for a pure superfluid used to obtain the frequency uj unita rity 
The decrease of the mode frequency which we predict as we approach the unitarity limit 
from the BEC side is thus qualitatively consistent with the experimental results in superfluid 
6 Li0,Q- 

The small damping of the breathing modes in superfluid 6 Li:2fj m &y also be related to 
the "suppression effect" we have discussed in this paper, i.e., collective modes always stay 
below the two-particle excitation threshold and thus are stable with no damping at T = 0. 
On the other hand, our results show that the frequencies of the collective modes smoothly 
change as we pass through the crossover region, and does not predict the observed sudden 
increase of the radial breathing mode frequency in the BCS regime of superfluid 6 Lij3|. A 
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breakdown of the superfluid phase or some other cause would seem to be the source of this 
sudden increase in damping, but this requires further study. 



VI. SUMMARY 



In this paper, using a collisionless HFB-GRPA type theory at T = 0, we have investigated 
the collective modes in a trapped superfluid Fermi gas with a Feshbach resonance. In 



essence, we have extended our previous work|8( on the collective modes in a uniform Fermi 
superfluid to include an isotropic harmonic trap. Our present calculation of the collective 
modes makes crucial use of our recent work|22| on the single-particle excitations of a trapped 
Fermi superfluid, given by the solutions of the Bogoliubov-de Gennes coupled equations. 
Using these eigenstates, we have calculated the gene ralized density correlation function in 
the HFB-GRPA formalism, which has been shown[32] to be consistent with both Goldstone's 
theorem and Kohn's theorem in the crossover region of a trapped Fermi superfluid. 

One of the major results of the present paper is that we have shown in detail how 
quadrupole and monopole mode frequencies are suppressed by the threshold energy of the 
two-particle continuum (twice the single-particle excitation gap E g ). When 2E g is smaller 
than the collective mode energy (= u hd ) predicted by a pure superfluid (or hydro dynamic) 
theory, the mode frequency can deviate significantly from uJh.d.- As we showed in Ref. |27[ 
(see also Fig. |3]of this paper), the single-particle excitation gap E g is of the order of the trap 
frequency Uq in a trapped Fermi superfluid in the crossover region. This is much smaller 
than the magnitude of the composite order parameter A(r = 0) at the center of the trap, 
which is comparable to the Fermi energy £p. This is because E„ is determined by the lowest 



energy Andreev bound state localized at the edge of the trap[27|, |28[. The suppression of 
the collective mode by the two-particle continuum should be observable, and would be even 
more important at temperatures close to T c . 

We have emphasized that the use of the T = equation of motion of a pure superfluid: 
(which assumes that the density of superfluid equals the total density) precludes any 
inclusion of effects of the two-particle continuum. Thus, the work presented here extends 
previous work. We also note that our work is built on the calculation (see Ref. Q and 
Appendix A) of the correct BdG BCS excitations of a trapped Fermi superfluid. This allow 
us to highlight the feature that the single-particle energy gap E g arises from low energy BCS 
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excitations, localized near the edge of the trap. These states, by definition, are left out any 
theory based on the local density approximation (LDA). 

In the BEC regime, where the Fermi chemical potential \x is negative, the two-particle 
excitation spectrum exhibits a large gap, with 2E g ~ 2\/j,\ ^> c^o- In this limit, the low-energy 
collective modes with uj ~ ujq are no longer much effected by the higher energy two-particle 
continuum. 

Our results in this paper may be useful as a way to obtain useful information about 
the single-particle excitation gap energy E g ~ u in the crossover region. Recently, the 
single-particle Fermi excitation spectrum was studied using the rf-tunneling spectroscopy 



in 6 Li 



24j . In this experiment, the energy gap is evaluated from an emergent peak in the 



spectrum as the gas becomes superfluid. This observed gap actually only gives information 
about the magnitude of the composite order parameter A(r = 0) in the center of the trap, 
which is of the order of the Fermi energy £f(^ > ^o) m the crossover region. In principle, 
such rf-tunneling experiments should be able to observe the true excitation gap at much 
smaller energies E g ~ ujq, but this would require very good energy resolution. Observing 
the suppression of a collective mode frequency discussed in this paper due to the presence 
of the two-particle continuum may be an alternative way to measure the magnitude of E g . 
Although we have assumed an isotropic trap in this paper, the same suppression will also 
occur in the anisotropic trap used in recent experiments. 

In this paper, we have only considered the case of T = 0. At finite temperatures, 
Bogoliubov quasi-particles are excited thermally, which will give rise to Landau dam ping of 



the collective modes. In addition, fluctuation effects from the Cooper-channel j 
become crucial, especially near T c . In future work, we hope to extend the present results at 
finite temperatures. 
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APPENDIX A: SINGLE-PARTICLE EXCITATIONS IN THE BCS-BEC 
CROSSOVER REGION 



In this appendix, we review some key results in Ref. 27[ for the single-particle excitations. 
We expand the fermion field operator ^(r) in the eigenfunctions fni m ( r ) = u ni( r )Yim(&) of 
the harmonic potential V t f ap (r), as 



(Al) 



nlm 

Here Yi m {9) is a spherical harmonic, and the radial wavefunction has the form 

u£(r) = v^ ) 3/4 y^ ^ 1/2) , e-^f'Li + V2 (f 2 ) (f s (A2) 

where L l + l l 2 {r 2 ) is the Laguerre polynomial. Substituting Eq. (jAl|) into -£^jfb m Eq. (|2.6|) . 
we obtain i^HFB = S/m-^HFB-; where the (/, m)-component H^ B is given by 

U 



-^HFB ~ ^ ^>nl C nlma C nlma n ^nn lC nlma C n'lma 

n,cr z nn',a 



^nn' [ C Lmt C n''-m| ^" C ' 



(A3) 



Here ^ = hco (2n + I + 3/2) — /i is the eigenenergy of /^ m (r). i^ n , and J^ n , represent 
the scattering matrix elements involving the composite pair potential A(r) and the Hartree 



potential — ^np(r), respectively. They are defined by 



PL* = / r 2 dru^r)A(r)u^ n (r), 
Jo 



V 



r 2 rfra^(r)n F (r)M^(r). 



(A4) 
(A5) 



The matrix element F l nn , in Eq. (jA4|) includes both intrashell pairing (n = n') and intershell 
pairing (n ^ n'). 

The quadratic Hamiltonian in Eq. (jA3|) can be diagonalized by the unitary transformation 



/ r \ 

Coiml 



c, 



t 

01 -ml 



w 



\ c Nil -ml ) 



^ lOlml X 



iNilm] 
lOlml 



\ iNilml ) 



(A6) 



21 



where the 2(N t + 1) x 2(N t + l)-orthogonal matrix W l is taken so that W l ^H^ B W l be 
diagonal. This requirement leads to the BdG coupled equations for the wavefunction W- n , 



where i, n can go from to 2(iV z + l) [see Eq. (2.25) of Ref. 27]]. In Eq. (jA6|) . 7j; mo - represents 
the annihilation operator of a Bogoliubov quasiparticle excitation, with eigenenergy E^. The 
resulting diagonalized H^ B has the form (constant terms are dropped) 

H^ FB = J2E^ jlmaljlma (*5>0). (A7) 

3=0 

As usual in the BCS theory, one needs a cutoff uj c = hu (N c + 3/2) to avoid ultraviolet 
divergences. The maximal radial quantum number Ni in Eq. (|A7|I is given by the largest 
integer bounded by (N c — /)/2. Using the solutions of the BdG equations, the Cooper-pair 



order parameter A 
the BdG equations 



r) and the Fermi atom density rip(r) are given by the solutions W- n of 



n 7 i 1 Ni 

A(r) = UY! ^—u nl {r)u nn {r) £ W l Rl+nA+j W l n , +lA+j . (A8) 

nn'l j=0 
2/ -|- 1 ^ l 

Mr) = E —TT- um(r)u n n(r) ]T [W^ +lA+j W^ +liNi+j + W l Rl+nd+1 W l Rl+nW ] . (A9) 

nn'l ^ n j=0 

Here iVj = Ni + 2 and the prime in J]' refers to the cutoff o; c . 

Since the molecular condensate component 0jw(r) is isotropic, it can be expanded in 
terms of the radial component «^(r)[57| of the wavefunction fn( m (r) = u^i{r)Yi m {9) for a 
Bose molecule in the harmonic potential V t ^[ p (r). We obtain 

</>m(t) = -iE a ^oW 

V 47T n 

= _gry <o(0 [°° dr'r ,2 A(r'), (AlO) 

where E% = hu (2n + I + 3/2) is the energy of the molecular state /^f m (r). The last 
expression in Eq. (|A10j) is obtained by using the identity in Eq. (|2.4jl . 

The boson mean- field Hamiltonian 5H M in Eq. (Q can be easily diagonalized, reducing 

to 

5H M = J2( E nl + ^ - ^Ijnlm, (All) 
nZm 

where b n i m is the annihilation operator of the molecular eigenstate f^ m (r). 
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To describe the BCS-BEC crossover region at T = 0, we extend the strong-coupling 
theory developed by Leggett jl^. to the case of superfluidity in a Fermi gas trapped in 
a harmonic potential with a Feshbach resonance. The key is that one needs to take into 
account the deviation of the chemical potential /i from the Fermi energy due to the formation 



of a bound state Bose condensate 



121 ]. This effect is incorporated into the theory at T = 



by solving the BdG equations together with the equation for the number of Fermi atoms. 
In the presence of a Feshbach resonance, a fraction of these atoms form molecular bosons. 
At T = 0, our discussion assumes all these molecules are Bose condensed, so that the total 
density of atoms is given by 



n(r) = 2|0 M (r)| 2 + n F (r) = 2n M (r) + n F (r). 



(A12) 



The equation of the state fi(N) giving fi as a function of the total number of atoms to 
determine \x is obtained by integrating Eq. (jA12|) over r. The result is 



Ar = 2 E«n + E(^ + 1 ) 



\W l - I 2 4- \W l - I 2 
I vv n+l,Ni+n'\ ' I vy Ni+n,n'+l\ 



(A13) 



nn'l 



The self-consistent solution W in and (|A13|) are the basis of our mean- field crossover theory 
at T = 0(22] . The resulting BCS excitations are used to calculate the response functions in 
the text of this paper. 



APPENDIX B: MOLECULAR BOSE PROPAGATOR 

The renormalized 2x2 matrix Bose Green's function describing molecules, which includes 
the Feshbach resonance effect within the GRPA, is given by|8( (we use a matrix notation for 
r and r') 



D(u) 



b°(u). 



(Bl) 



The explicit form for D°(u) is given in Eq. 1)3.6)1 . The self-energy of the molecules S(a>) is 
due to the fluctuations associated with the particle-particle fluctuations in the Fermi gas, 



2 ' 



irv 



p. 



(B2) 



Here we have introduced a projection operator rj[A] which extracts the (1,1), (1,2), (2,1), 
and (2,2) components from a 3 x 3-matrix A. The 2x2 matrix P is defined in Eq. (13.5)1 . 
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Excitations from the molecular condensate are determined from the poles of the renor- 
malized Bose Green's function D(uj) in Eq. (jBlj) . This has been discussed in Ref.jsJ for the 
case of a uniform superfluid Fermi gas. The excitation energy of the dressed molecule is 
obtained from the zeros of the determinant 







det 
det 
det 



1 _ ^no^pifj 

i - f BMnV 
u 



i + -nV) 



p 



i + |ftv; 



i h -nV) - f ^MnV) 



2 v 



det 


1 + | 


U-g 2 t B{u) 


n°(o;) 


det 


1 + f nv) 





(B3) 



Here 1 + |[Z7 — ^-B(u;)]ri (u;) and 1 + ^fl (u;) in the last line are 3 x 3-matrices in the 
(pi) P2, Pz) -space introduced in Sec. II. Comparing Eq. ()B3j) with Eq. ()3.3|) . we see that 
Bose molecular excitations are the same as the collective mode resonances of the (generalized) 
density correlation functions [unless the denominator in Eq. (jB3|) diverges] . This equivalence 
is to be expected since the Feshbach coupling in the Hamiltonian in Eq. ()2.1|) mixes the 
molecular and Cooper-pair channels. 
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FIG. 1: Single particle density of states N(ui) in the weak-coupling BCS region (where the chemical 
potential /j, is positive). The rapid oscillations in the spectrum are the discrete Bogoliubov quasi- 
particle energies in a harmonic potential. In this figure, we have introduced a small imaginary part 
T = 0.005cjo to slightly broaden the Bogoliubov quasi-particle peaks. See also Ref. [2^ . 

FIG. 2: Spatial variation of the composite order parameter A(r) in the BCS region. The pa- 
rameters are the same as Fig. ^ The dashed line shows the "diagonal" potential measured from 
the chemical potential fi, given by Eq. (|4.4j) . The quasiparticles also feel the Hartree potential 
-|n(r), but this is not included. 



FIG. 3: (a) Magnitude of the composite order parameter A(r = 0) at the center of the trap 
and the single-particle excitation gap E g in the BCS-BEC crossover region, (b) Fermi chemical 
potential fi in the crossover region. The dotted line in panel (a) shows the absolute value of the 
chemical potential \fi\ when [i < 0. 
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FIG. 4: (a) Local density of states N(u, r) defined by Eq. (H2J) in the BCS region {v/ey = 2). In 
this case, A(r = 0) = Wloq S> c<jo at the center of the trap. We take U = 0.52£F) with the other 
parameters the same as in Figs. ^andEl To smooth the sharp delta-functions, we have introduced 
a small imaginary part 7 = 0.2ujq in the Bogoliubov excitation energies, (b) N(u, r) based on the 
LDA. In this figure, we take 7 = O.lwo- The ridge in the LDA single-particle spectrum occurs at 
uj = A(r). 

FIG. 5: (a) Local density of states N(u,r) in the BCS region (v/ey = 2) in the case where 
A(r = 0) = 1.86u;o at the center of the trap. We take U = 0.35£f- As usual, the finite magnitude 
of N(uj = 0, r) is due to the small imaginary part 7 = 0.2loq introduced into the Bogoliubov quasi- 
particle excitations, (b) N(uj,r) in the strong-coupling BEC region, where the Fermi chemical 
potential \i is negative. 

FIG. 6: Frequency of the quadrupole mode (L = 2 and n = 0) in the BCS-BEC crossover, for 
two different values of v. The threshold energy 2E g of two-particle continuum is evaluated directly 
from the calculated energy gap in N(uj), as shown in Fig. 1. 

FIG. 7: Effective interaction between molecules mediated by free Fermi atoms in the BEC limit. 
The solid and dashed lines represent the Fermi and molecular Bose Green's functions, respectively. 

FIG. 8: Spectrum of the strength function 633(0;, L = 2) for U = 0.35£f- The lowest peak 
corresponds to the quadrupole mode. The arrows show where the two-particle spectrum (u > 2E g ) 
starts. In this figure, for clarity, we have given the Bogoliubov energies a small imaginary part 
7 = 0.005u;o. The same peak broadening is also used in Figs. andlT2l 

FIG. 9: The same as in Fig. El for U = 0.52£f- Since the two-particle threshold 2E g ~ 2.5wo [see 
Fig. 6(a)] is much higher than the trap frequency loq in the BCS and crossover region, only the 
quadrupole mode peak appears in the low energy spectrum. 

FIG. 10: Frequency of the monopole mode (L = and n = 1) in the BCS-BEC crossover. 2E g is 
the threshold energy of the two-particle continuum. 
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FIG. 11: Spectrum of the strength function 633(0;, L = 0) for U = 0.35£f- The lowest energy 
peak corresponds to the monopole mode. Since the spectral weight of the monopole is very weak 
in panel (a), an inset shows a magnified spectrum around this peak. In panel (a), the strong peak 
at <jj ~ u>q comes from two-particle excitations. The arrow shows the two-particle threshold energy 
2E g . The expected singularity around to = arising from the Goldstone mode is not shown in this 
figure. 

FIG. 12: The same as Fig. HUfor U = 0.52e F . 
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